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Abstract 

In the present work we study non-thermal leptogenesis and baryon 
asymmetry in the universe in different neutrino mass models discussed 
recently. For each model we obtain a formula relating the reheating 
temperature after inflation to the inflaton mass. It is shown that 
all but four cases are excluded and that in the cases which survive 
the inflaton mass and the reheating temperature after inflation are 
bounded from below and from above. 



1 Introduction 



The Standard Model (SM) of particle physics dj is a very successful theoreti- 
cal framework for all low-energy phenomena. However, it is widely considered 
to be a low-energy limit of some underline fundamental theory. Perhaps the 
most direct evidence for physics beyond the SM is the recent discovery that 
neutrinos have small but finite masses [21 El H] • A simple and natural way to 
explain the tiny neutrino masses is via the seesaw mechanism [5 J . According 
to that, the existence of super-heavy right-handed neutrinos is postulated 
and the smallness of the masses of the usual SM neutrinos is due to the 
largeness of the masses of the new neutrinos. Solar, atmospheric, reactor 
and accelerator neutrino experiments (for a summary of three-flavour neu- 
trino oscillation parameters see e.g. seem to indicate neutrino masses 
in the sub-eV range (0.001 eV < m v < 0.1 eV), which implies that heavy 
right-handed neutrinos weigh ~ 10 10 GeV — 10 15 GeV [7j. 

On the other hand, the baryon asymmetry in the universe (BAU) is one of 
the most challenging problems for modern cosmology. Both Big-Bang Nucle- 
osynthesis jH] and CMB data (for example from WMAP 0) show that in the 
universe one baryon corresponds approximately to one billion photons. This 
very small number should be computable in the framework of the theory of 
the elementary particles and their interactions we know today. Nowadays, the 
most popular way to obtain the BAU is through leptogenesis ^U]- Initially 
a lepton asymmetry is generated through the out-of-equilibrium decays of 
right-handed neutrinos and then the lepton asymmetry is partially converted 
to baryon asymmetry through the non-perturbative "sphaleron" effects [TT] . 
In general leptogenesis can be thermal or non-thermal. Thermal leptogenesis 
usually requires very high reheating temperature after inflation ^2]- This 
can be problematic because of the gravitino constraint. In supersymmetric 
models (for reviews in supersymmetry see e.g. ^3] and for supersymmetry 
in cosmology see e.g. [14j ) with spontaneous supersymmetry breaking the 
superpartner of the graviton, the gravitino, gets a mass depending on how 
the supersymmetry is broken. In gravity mediated supersymmetry break- 
ing the gravitino mass is in the range m 3 / 2 — 100 GeV — 1 TeV and the 
gravitino (if not the lightest supersymmetris particle) is unstable with a life- 
time larger than Nucleosynthesis time ~ 1 sec and dangerous for cosmol- 
ogy. This gravitino problem can be avoided provided that the reheating 
temperature after inflation is bounded from above in a certain way, namely 
T R < (10 6 - 10 7 ) GeV TO]. 
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Therefore one can see that heavy right-handed neutrinos can have impor- 
tant implications both for particle physics and cosmology. Various neutrino 
mass models ^3 QB] have been proposed and their predictions on neutrino 
masses and mixings have been studied thoroughly. The requirement for the 
right baryon asymmetry in the universe as well as for the right phenomenol- 
ogy for light neutrino masses and mixings puts severe constraints on right- 
handed neutrinos. Recently six concrete neutrino mass models were discussed 
and a comparison of numerical predictions on baryon asymmetry for these 
models was presented [TU] . Two of the models were almost consistent with 
the observed BAU, while the rest of them predicted either a small (rj < 10~ 19 ) 
or a large (r] > 10 -6 ) baryon asymmetry. The analysis was performed in the 
framework of thermal leptogenesis. The aim of the present work is to study 
the same models in the framework of non-thermal leptogenesis and derive the 
constaints on the inflaton mass and the reheating temperature after inflation. 

Our work is organized as follows. After this introduction we review the six 
neutrino mass models in section 2 and we discuss non-thermal leptogenesis 
for these models in the third section. Our results are presented in section 4 
and we conclude in the last section. 

2 Review of the different neutrino mass mod- 
els 

Here we give a brief review of the six neutrino mass models ^H] discussed 
recently in [19J. The interested reader can find more details in ^HJ HH]- In 
particular, all the information about the models are collected in Appendix 
A of [T^j. There is one normal hierarchical model (NHT3), two inverted 
hierchical models (InvT2A, InvT2B) and three degenerate models (DegTlA, 
DegTIB, DegTIC). According to seesaw mechanism, the light left-handed 
neutrino mass matrix m u , the heavy right-handed neutrino mass matrix Mr 
and the Dirac neutrino mass matrix are related as follows 

m v = m D M R ^ , mF D (1) 

where M^ 1 is the inverse of Mr and is the transpose of mo- The pre- 
dicted values of the neutrino mass-squared differences and mixing parameters 
are shown in table 1. 
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Type 


Am 2 21 [10- 5 eV 2 


Aml 3 [10- 3 eV 2 


tan 2 9 12 


sin 2 26> 23 


sin 6»i3 


DegTIA 


8.80 


2.83 


0.98 


1.0 


0.0 


DegTIB 


7.91 


2.50 


0.27 


1.0 


0.0 


DegTIC 


7.91 


2.50 


0.27 


1.0 


0.0 


InvT2A 


8.36 


2.50 


0.44 


1.0 


0.0 


InvT2B 


9.30 


2.50 


0.98 


1.0 


0.0 


NHT3 


9.04 


3.01 


0.55 


0.98 


0.074 



Table 1: Predicted values of the solar and atmospheric neutrino mass-squared differences 
and three mixing parameters (from 



Type 


Case (i): \Mj\ 


Case (ii): \Mj\ 


DegTIA 


4.28 x 10 9 , 1.16 x 10 10 , 


3.47 x 10 7 , 9.42 x 10 7 , 




3.84 x 10 13 


3.81 x 10 13 


DegTIB 


4.05 x 10 7 , 6.16 x 10 11 , 


3.28 x 10 5 , 4.98 x 10 9 , 




7.6 x 10 13 


7.6 x 10 13 


DegTIC 


4.05 x 10 7 , 6.69 x 10 12 , 


3.28 x 10 5 , 4.85 x 10 11 , 




6.99 x 10 12 


7.81 x 10 11 


InvT2A 


3.28 x 10 8 , 9.70 x 10 12 , 


2.64 x 10 6 , 7.92 x 10 10 , 




6.79 x 10 16 


6.70 x 10 16 


InvT2B 


5.6527 x 10 10 , 5.6532 x 10 10 , 


4.5971 x 10 8 ,4.5974 x 10 8 , 




5.38 x 10 16 


5.34 x 10 16 


NHT3 


6.51 x 10 10 , 7.97 x 10 11 , 


5.27 x 10 8 ,6.45 x 10 9 , 




1.01 x 10 15 


1.01 x 10 15 



Table 2: The three right-handed Majorana neutrino masses in GeV (from |19|). 



Type 


Case (i): 


e 


Case (ii): 


e 


Case (i): 


V 


Case (ii): 


V 


DegTIA 


2.10 x 10 


-6 


1.71 x 10 


-8 


4.99 x 10 


-9 


4.06 x 10" 


-ii 


DegTIB 


2.66 x 10' 


-18 


2.16 x 10- 


-20 


1.60 x 10- 


-23 


1.30 x 10- 


-25 


DegTIC 


1.74 x 10" 


-18 


1.69 x 10- 


-20 


1.05 x 10" 


-23 


1.02 x 10- 


-25 


InvT2A 


1.59 x 10" 


-14 


1.27 x 10- 


-16 


9.94 x 10" 


-19 


7.96 x 10- 


-21 


InvT2B 


1.47 x 10 


-2 


1.62 x 10 


-4 


5.40 x 10 


-5 


5.94 x 10 


-7 


NHT3 


5.90 x 10 


-7 


4.78 x 10 


-9 


2.17 x 10 


-9 


1.76 x 10- 


-11 



Table 3: Calculation of CP asymmetry e and baryon asymmetry 77 for each neutrino 
mass model (from 
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In thermal leptogenesis for the SM case the BAU r/ = n^/n 7 = 6.1 x 10 10 is 
computed by the formula ^Hl 

tj = 0.0216«;e (2) 
where the CP asymmetry e is defined as 

r - f 

e= r+T < 3 > 

with r = r(A r i — * li,(jy) and V = T(N\ — > lL<j>) the decay rates, while the dilution 
factor k is determined by numerical integration of Boltzmann equations. However 
it can be estimated by ^H] 

(4) 



for < K < 10 and by 



2V9 + K 2 
0.3 



'" K(lnK) oe ^ 

for 10 < K < 10 6 , with K the decay parameter K = rhi/m*, where m* is the 
equilibrium neutrino mass m* = 1.08 x 10 -3 eV and rh\ is the effective neutrino 
mass defined as 

with v the electroweak scale, M\ the mass of N\ and h the matrix for the neutrino 
Yukawa couplings. The three right-handed neutrino masses for each model are 
shown in table 2 while the CP asymmetry and baryon asymmetry are shown in 
table 3. The Dirac neutrino mass matrix mo can be either the charged lepton 
mass matrix m; (case (i)) or the up-quark mass matrix m u (case (h)). We see that 
NHT3 and DegTlA models are almost consistent with the observed BAU, while 
the rest of the models lead either to very small baryon asymmetry, rj < 10~ 19 
(DegTIB, DegTIC, InvT2A), or to large baryon asymmetry, rj > 10~ 6 (InvT2B). 



3 Non-thermal leptogenesis 

We start by introducing three heavy right-handed neutrinos (one for each family) 
Ni, i = 1, 2, 3 with masses Mi, M2, M3, which interact only with leptons and Higgs 
through Yukawa couplings. In supersymmetric models the superpotential that 
describes their interactions with leptons and Higgs is [2111 

W l = Y ia N % L a H u (7) 

where Yi a is the matrix for the Yukawa couplings, H u is the superfield of the 
Higgs doublet that couples to up-type quarks and L a (a = e, [i, r) is the superfield 
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of the lepton doublets. Furthermore, we assume that after the slow-roll phase 
of inflation, the inflaton decays dominantly to right-handed neutrinos through 
Yukawa couplings and for supersymmetric models the interaction is described by 
the superpotential [2J 

W 2 = ]T \ t SNf Ni (8) 

i 

where Aj are the couplings for this type of interaction and S is a gauge singlet 
chiral superfield for the inflaton. With such a superpotential the inflaton decay 
rate is given by |U 

r (j} = r(cj ) ^N l N l ) = ^-\\ i \ 2 M I (9) 



The reheating temperature after inflation Tr is given by [22] 

(i> M vl fi 2 (io) 



\47r 3 5r !f! 

where M p i is Planck mass and g* is the effective number of relativistic degrees of 
freedom at the reheating temperature. For the reheating temperatures that we 
shall consider all the particles are relativistic and for MSSM = 915/4 = 228.75, 
while for SM g* = 427/4 = 106.75. 

Any lepton asymmetry Yl = tll/s produced before the electroweak phase 
transition is partially converted into a baryon asymmetry Yb = n B /s via sphaleron 
effects [TT] • The resulting Yb is 

Y B = CY L (11) 

with the fraction C computed to be C = -8/15 in the MSSM and C = -28/79 
in the SM [2H1- The lepton asymmetry, in turn, is generated by the CP- violating 
out-of-equilibrium decays of the heavy neutrinos 

N^IH*, N -> IH U (12) 



In the framework of non-thermal leptogenesis the lepton asymmetry is given by [24 

Y L = ^BR(<t>^N 1 N 1 )^e (13) 

where Mj is the inflaton mass, Tr the reheating temperature after inflation, e the 
CP asymmetry and BR is the branching ratio for the decay of the inflaton to the 
lightest heavy right-handed neutrino. The decay is kinematically allowed provided 
that 

Mi > 2Mi (14) 
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We will assume that BR ~ 1, that is the inflaton decays practically only to the 
lightest of the right-handed neutrinos. This is possible even if the inflaton is 
heavy enough to decay to all right-handed neutrinos as long as |Ai| 2 3> | A2 1 2 , | A3 1 2 . 
Combining the above formulae we obtain 

or 

t »={W m ' (16) 

From the WMAP data jU] we know that 

rj B = M = 6.1 x 10- 10 (17) 

If we recall that the entropy density for relativistic degrees of freedom is s = 
h e jf^-T 3 and that the number density for photons is n 7 = 2< ^P T 3 , one easily 
obtains for today that s = 7.04n 7 . Thus for Yb we have 

Y B = 8.7 x KT 11 (18) 

Following [23j we shall consider that M\ > IOOTr, because in that case the neutrino 
N\ is always out of thermal equilibrium. Finally we recall that Mj > 2M\. 



4 Results 

Now we can present our results. We shall begin with the SM case first and we shall 
use for the fraction C the SM value, namely C = —28/79. For each neutrino model 
(12 cases in total) the CP asymmetry e as well as the right-handed neutrino mass 
Mi are known. Therefore we have i) a formula relating the reheating temperature 
to the inflaton mass, ii) a lower bound for the inflaton mass Mj > 2M%, and iii) 
an upper bound for the reheating temperature Tr < 0.01 M\. Furthermore, using 
the relationship between Tr and Mi we are able to convert the upper limit for 
Tr to a corresponding upper limit for Mj and also the lower limit for M\ to a 
corresponding lower limit for Tr. So both Tr and Mi are bounded both from 
above and from below. Let T™ n and T^ ax be the lower and higher value for the 
reheating temperature respectively. Then T™" < Tr < T^ ax and obviously it is 
required that T^ ax > T™ n , which is not satisfied for all cases. In fact most of the 
cases are excluded. The only cases for which the constraint is satisfied are: 
- DegTIA, case (i), for which: 

8.56 x 10 9 GeV < Mi < 5.49 x 10 11 GeV (19) 
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6.67 x 10 5 GeV <T R < 4.28 x 10 7 GeV (20) 
- NHT3, case (i), for which: 

1.3 x 10 11 GeV < Mi < 2.35 x 10 12 GeV (21) 

3.6 x 10 7 GeV < T R < 6.51 x 10 8 GeV (22) 
-InvT2B, case (i), for which: 

1.13 x 10 11 GeV < Mj < 5.09 x 10 16 GeV (23) 

1.25 x 10 3 GeV < T R < 5.65 x 10 8 GeV (24) 
-InvT2B, case (ii), for which: 

9.2 x 10 s GeV < Mj < 4.55 x 10 12 GeV (25) 

9.29 x 10 2 GeV < T R < 4.6 x 10 6 GeV (26) 

One can see from the results presented above that inflationary models in which 
Mj ~ 10 13 GeV, like e.g. chaotic [25] or natural |26j inflation, are compatible 
only with one neutrino model (InvT2B, case (i)). Furthermore, for a concrete 
inflationary model with a given inflaton mass our results allow us to know what 
the reheating temperature must be and also what the inflaton decay rate is 
and what the inflaton Yukawa coupling |Ai| is. For example, in chaotic or natural 
inflation we obtain 

Mj ~ 10 13 GeV (27) 
T R ~ 10 5 GeV (28) 
T ~ 10~ 8 GeV (29) 
| Ai| ~ 10~ 10 (30) 

At this point we should add a comment regarding the gravitino constraint. If we 
add supersymmetry in order to address the gravitino problem, then the expression 
for the baryon asymmetry in MSSM will change slightly by a numerical factor of 
order one. So we could use the results obtained so far for the SM case. If we require 
that T R < (10 6 - 10 7 ) GeV then we see that the models InvT2B and DegTIA are 
already compatible with the gravitino constraint, the model NHT3 is marginally 
compatible (for the lower values for T R ) with the gravitino constraint and finally 
the model InvT2B can be made compatible with the gravitino constraint lowering 
the upper bound for T R 

1.25 x 10 3 GeV <T R < (10 6 - 10 7 ) GeV (31) 
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5 Conclusions 



In the present work we have studied non-thermal leptogenesis in six neutrino mass 
models proposed earlier and discussed recently in the literature. For each model 
we have obtained a formula relating the inflaton mass Mj to the reheating tem- 
perature after inflation Tr. In fact according to this formula Tr is proportional 
to Mj. Hence, the bigger the inflaton mass the bigger the reheating tempera- 
ture. In a concrete inflationary model (chaotic |25| . natural 26 , supersymmetric 
hybrid [2Zj etc) with a given mass for the inflaton, the right baryon asymmetry 
implies a certain reheating temperatute after inflation. This in turn implies a cer- 
tain decay rate for the inflaton field and a certain value for the inflaton Yukawa 
coupling. Furthermore, kinematical reasons and the requirement for non-thermal 
leptogenesis lead to a lower and an upper bound both for Mi and Tr. Our results 
show that in most of the neutrino models under study the lower bound is not 
compatible with the upper bound and therefore only four cases survive. If we also 
take into account the gravitino constraint Tr < (10 6 — 10 7 ) GeV, then in one of 
these cases the reheating temperature is even more constrained. 
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